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Various algebraic structures of degenerate four-wave mixing equations of optical phase conjuga- 
tion are analyzed. Two approaches (the spinorial and the Lax-pair based), complementary to each 
other, are utilized for a systematic derivation of conserved quantities. Symmetry groups of both 
the equations and the conserved quantities are determined, and the corresponding generators are 
written down explicitly. Relation between these two symmetry groups is found. Conserved quanti- 
ties enable the introduction of new methods for integration of the equations in the cases when the 
coupling r is either purely real or purely imaginary. These methods allow for both geometries of the 
process, namely the transmission and the reflection, to be treated on an equal basis. One approach 
to introduction of Hamiltonian and Lagrangian structures for the 4WM systems is explored, and the 
obstacles in successful implementation of that programe are identifled. In case of real coupling these 
obstacles are removable, and full Hamiltonian and Lagrangian formulations of the initial system are 
possible. 



I. INTRODUCTION 



There is a short story prefacing the paper. The work 
on symmetries in optical phase conjugation started in 
early nineties by Predrag Stojkov and Milivoj Belie. It 
was interrupted by Predrag's leaving for America in 1992. 
During the stay of M. Belie at the Texas A&M Univer- 
sity in 1995 and 1996, the problem and an early draft of 
the paper were brought to Marko's attention. At the time 
he was phasing out of quasicrystals, and was open to new 
ideas. Marko liked the problem and agreed to participate. 
He read the manuscript, made numerous changes, and 
suggested a new direction to it. Owing to his commit- 
ments for the sabbatical at Cornell University in 1995 
and the visit to Israel in 1996, it was decided to postpone 
the serious work after he is back. However, during the 
Israeli visit Marko was diagnosed with the brain tumor. 
The paper is left essentially unchanged. It is dedicated to 
his memory. 

Steady-state four-wave mixing (4WM) equations de- 
scribing optical phase conjugation (OPC) in photorefrac- 
tive (PR) crystals have been solved up to now in a num- 
ber of ways . A common feature of all solution meth- 
ods is that, first, conserved quantities are determined, 
and then the number of equations is reduced. However, 
the determination of conserved quantities and the reduc- 
tion of equations is usually performed in an ad hoc man- 
ner. Furthermore, the solution of the OPC equations in 
the two basic geometries of the process, the transmission 
geometry and the reflection geometry, is usually obtained 
using unrelated methods. 

Apparent symmetries of wave equations have not been 
used up to now to facilitate the analysis and the 



solution of the problem. In this work the symmetries 
of the equations and the integrals of motion are investi- 
gated and used to present a unified method for systematic 
derivation of conserved quantities, an equal treatment of 
both geometries, and the reduction in the number of in- 
dependent variables Such an analysis allows for not 
only an easier handling of otherwise cumbersome and un- 
related relations, but also for a deeper understanding of 
the physics of the process. Also, rudiments of a formal 
presentation of the problem along the lines of the theory 
of dynamical systems are presented. 

The geometry of the process is simple. Three laser 
beams intersect within a piece of the PR crystal: two 
counterpropagating laser pumps Ai and A2, and a signal 
A4. Owing to the PR effect, a fourth wave A^ is gener- 
ated inside the crystal, that counterpropagates to, and is 
the phase conjugate replica of the signal. There are two 
main channels along which the generation may proceed. 
In the first one, the signal wave builds a diffraction grat- 
ing with the pump Ai. The other pump is diffracted off 
that grating and transmitted across the crystal into the 
PC wave A3. This is the so-called transmission geometry 
(TG) of the process. 

In the second channel the signal interferes with the 
pump A2, and the beam Ai is reflected off the grating 
into the PC wave. This is the reflection geometry (RG) 
of the 4WM process. It is assumed that all waves oscil- 
late at the same frequency (the degenerate 4WM). Also, a 
steady-state is assumed, and all beams are approximated 
by plane waves. 

The equations of interest are the slowly varying enve- 
lope wave equations describing 4WM in PR media |Q] . In 
TG, they are of the form 

IA[ = VQtA^, 



1 



lA' 



-TQtAi, 



(1.1) 



where / — l^iP the total intensity, F is the 

coupling constant (complex in general, but often real in 
PR media), Qt — A1A4 + A2A3 represents the diffrac- 
tion grating amplitude for TG, the prime denotes spatial 
derivative along the propagation z direction, and the bar 
denotes complex conjugation. 

In RG, the equations are given by 



IA[ = 


-TQrAs, 




IA'2 = 


-TQRAi, 




M'3 = 


-TQrAi, 




lA'^ = 


-rQRA2, 


(1.2) 



where the RG grating amplitude is given by Qr = 
A1A3 + A2A4. 

In this paper both geometries are treated on an equal 
footing, using a unified RG-like notation: 



Bi 
B2 

B3 
Bi 



A2, 

A^n^ + A^n^^ 

A^n^ + AiU^^ 



Ai in TG, 
A3 in RG, 
A3 in TG, 
A4 in RG, 



(1.3) 



where a is the switching variable, that has value +1 for 
TG and —1 for RG, and Il±a = (1 i cr)/2 are the corre- 
sponding "projectors." 

The "equations of motion" (EOM) are now 



IB[ = 
IB'2 = 

ibI = 

IB'4 = 



arQB3, 
aTQBi, 

-TQB2, 



(1.4) 

_B, P and the 



where the intensity is given by / = X]i=i 
grating amphtude hy Q — B1B3 + B2B4. 

The analysis is organized as follows. Two methods 
to derive the integrals of motion (lOM) are discussed in 
Section ||. The first method is based on the observation 
that 4WM EOM have a special symmetric form that is 
allowing an equivalent spinorial formulation. Such a form 
of EOM leads directly to the derivation of the full set of 
"regular" lOM as suitable bi-spinorial combinations. The 
symmetries of these lOM are the special unitary groups: 
SU{2) for TG, and SU{1,1) for RG. Initial spinor-like 
doublets of fields turn out to transform as the funda- 
mental irreducible representations of these groups, thus 
justifying the name "spinors". It is indicated how they 
can be used to reduce the number of dynamical variables. 

In the second method the Lax pair approach is uti- 
lized. The diadic products of the 4WM spinors are used 
as possible choices for the evolving member (£) of the 
Lax pair problem. The traces of products of these ma- 
trices represent lOM. It was established that all higher 



order lOM are various combinations of the basic lOM 
already obtained by the spinorial approach. At the end 
of Section || two special cases (F S M and F G iM.) are 
considered in some detail. 

In Section |l| attention is focused on the derivation 
of the symmetry groups of EOM, and the corresponding 
generators. The relation between these symmetries and 
the symmetries of lOM is discussed. 



In Section IV the symmetries of lOM are used to write 
the solutions of EOM (for F € R) in terms of elementary 
transcendental functions. Then an alternative solution 



procedure is explored. In the last part of Section [V the 
F g zM case is solved completely. 

The possibility of introducing the Hamiltonian and La- 
grangian description of 4WM EOM is explored in Section 
0. Section VI offers some conclusions and identifies open 



questions for future research. 



II. INTEGRALS OF MOTION AND THEIR 
SYMMETRIES 

A. Preliminaries 

In this work the 4WM equations are treated as a dy- 
namical system defined on the phase space V := V/p = 
E^, with the time variable z. Here V is the full sixteen- 
dimensional space ^ := = with the complex co- 
ordinates {x^} = {Bi,Bi} and p is the equivalence rela- 
tion (analyticity condition) satisfied by the 4WM system: 
^ (^^iy p ^ (^^)*] foj. ^ ^ xhe tangent 

space TV (the space of vector fields on V) , is spanned by 
the coordinate basis {9^} = {di = d/dBi\di = d/dBi}. 
The dynamics on the space V is given by the trajectory 
c : 'Rz ^ V . It is described by the velocity vector-field 
Fe TV: 

F = p{<jB3di - 3284) + p{aB4d2 - Bid3) + c.c. 

(2.1) 

where p := TQ/I. For a general function f{z,x) E 
C^(Mz X V), the corresponding evolution equation is 



(a. + F)f. 



(2.2) 



In general, an integral of motion (lOM) q{z, x) is a 
function that is constant along the trajectory c, i.e. 
{dz + F)q\c = {on-shell constancy). Here the more 
restrictive definition of lOM is used: instead of an 
on-shell constancy, the condition of off-shell constancy 
[{dz + F)q = in whole V] is used. Also, only the inte- 
grals q{x) without the explicit time-dependence are con- 
sidered, leading to the defining equation 



Fq = 0. 



(2.3) 
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There is no general procedure for finding lOM, and one 
has to take into account various specifics of the system at 
hand. For example, one may resort to the brute-force so- 
lution procedure, which is based on the observation that 

is a linear differential operator of the zeroth degree of 
homogeneity (i.e. upon its action on some homogeneous 
function Ps{x) of degree s, it produces another homoge- 
neous function of the same degree s) . This allows one to 
replace Eq. (|2.3|) by the set of infinitely many equations: 



Fqis) = 0, 



(2.4) 



for s S N, where are the components of q with the 
fixed degre e of homogeneity s {q = J2^i1is))- gen- 
eral, Eqs. {2A) can be solved by a general ansatz (the 
summation over the repeated indices is assumed): 



9(1) 
9(2) 



(2.5) 



which turns Eq. (2.4) into a set of conditions for the ma- 
trices a. After some algebra, one finds that there are no 
integrals of the first degree, and that there are several of 
the second and higher degrees. 



B. Spinorial formalism 

There exists a more elegant way to find integrals of 
motion of the second degree of homogeneity in the 4WM 
equations |^ . It is based on the fact that convenient pairs 
of columns ("spinors") can be formed: 



IV-i) = 



Now the equations of motion (|l.4| ) can be written as a 
pair of matrix equations: 

IV-.y^mlV-,), (j-1,2) 
where the "evolution matrix" m is 

(T/i 

-fl 



(2.6) 



(2.7) 



(2.8) 



and fi = TQ /I. The matrix m is traceless and Hermitian 
(for TG) or skew-Hermitian (for RG), so it belongs to the 
su{2) algebra for TG, or to the su{l, 1) for RG. 
lOM arc found using a simple Lemma: 

• Lemma 1 .• A pair of linear matrix equations 

= niyi IV'a) , IiAb)' = , (2.9) 

has an integral of motion {-tp a\ li' b) , if there exists a 
constant matrix n such that 



nms +m' A 



(2.10) 



where the dagger denotes the adjoint matrix. 

In 4WM the TOM are searched for in two possible 
forms, as or as (V^ij n | ijjj ) . For the first form, 

we have the defining equation (2.10) specified as 



nm + m^n = 0, 



(2.11) 



whereas for the second form of integrals, the defining 
equation is 



T 

m n 



0. 



(2.12) 



For the general (complex) F, the unique solutions (up to 
rescaling by a constant) of these defining equations are 
the matrices 











ni = 




I) 









for (2.11) and 



n2 = 



1 for TG, 
CT3 for RG, 



= -cna2 = 



(2.13) 



-ia2 for TG, 
1(72 for RG, 

(2.14) 



for ( 2.12 ), where aj are the Pauli matrices. The corre- 
sponding integrals are 



qi = {i)i\ni\^i) = h + ah, 

q2 = {ip2\ni\ip2) = h +(Th, 

qs = (V^2| ni \ipi) = BiBi - B3B2, 

q4 = (V'll "2 IV'2) = B1B2 + (JB3B4. 



(2.15) 



Since these lOM are present for arbitrary complex cou- 
pling F, they are said to be the regular lOM of the 4WM 
system. Later it will be shown that for special choices of 
F this system possesses additional (exceptional) lOM. 

Not all of the conserved quantities gi , 52 , 93 , 93 , 94 
and 94 are independent. There exists a relation 



1941 



cr|93| 



91 92 



(2.16) 



that reduces the number of (real) integrals of motion to 
five. Using the integrals, one can express the conjugated 
fields as dependent variables: 



Bi = [aBsqs + B2qi]/q4, 

B2 = [Biq2 - aB4qz]/q4, 

B3 = [^491 - ^l93]/94, 

B4 = [S293 + S392]/94- 



(2.17) 



A more natural way to reduce the number of variables 
using conserved quantities is to introduce polar coordi- 
nates, suggested by the form of the conserved quantities 
qi and 92: 



B4 = 
B2 = 
B3 = 



q^c{cr,ai) exp(i^i), 
q^cia, 02) exp(i/?2), 
q^s{a,ai)exp{iji), 
gi's(cr,a2)exp(i72). 



(2.18) 



3 



Here the new variables are the six angles: ai, a2, Pi, /32, 
7i, and 72. The symbols c and s stand for the trigono- 
metric cosine and sine functions in the TG case, and for 
the hyperbolic cosine and sine in RG case (see Appendix 
A). 

To further ilucidate the connection between lOM and 
the new variables, it is convenient to employ the integral 



-<j{i/jf,\m, 



(2.24) 



q, introduced below [Eq. (2.41)]. It can be expressed in 
terms of the new variables, and the result is 



9? 



92 + 29ig2[c(cr, 2q;i)c(o-, 2a2)+ 
-as{(7, 2ai)s(cr, 2^2) cos($)], (2.19) 



where $ = /3i+/32 — 7i~72 is the so-called relative phase. 
From the spherical and the hyperbolic trigonometry 
it is known that the expression in brackets can be un- 
derstood as a cosine (hyperbolic cosine) of some angle p, 
so that 2ai , 2q!2 and p are the sides of a spherical (hy- 
perbolic) triangle, and $ is its central angle. Therefore 



c(cr, p) 



q-qt 



ql 



const. 



2qiq2 

and $ only depends on ai and a2 '■ 

c{cr,p) - c{(T,2ai)c{a,2a2) 

cos(3>) = — - — — — - — 7 . 

as{<7, 2ai)s((T, 2a2) 

It is easy to check that: 

2 



(2.20) 



(2.21) 



|94| 



qiq2c 



ksP = qms (a, ^ 



(2.22) 



in agreement with Eq. (2.16). Thus, there are five in- 
dependent real conserved quantities: qi, 52, P, and the 
phases of qs and 54. The solution of EOM using these 
quantities is performed in Sec. IV. 



C. Lax Pairs 

The Lax pair representation (if it exists) helps deter- 
mination of the integrals of motion. In general, if given 
dynamical system admits a Lax pair representation 

^^[M,C], (2.23) 

where C and A4 are suitably chosen operators or matri- 
ces, then all the traces Tr{C'') {k € N) are lOM. The 
determination of such a Lax pair of operators {M , C) is 
usually the hardest part of the problem. The brackets in 
Eq. (2.23) stand for the commutator. 

In the case of 4WM, the suitable matrices are easy 
to find, starting from the compact form of the spinorial 



EOM (2.7) and their conjugated equations: 



where the index p — i,i and ji/'i) := nin2|'!/'i). The fol- 
lowing matrices 



satisfy the Lax pair equations 

The corresponding Laxian lOM are given by 



9mi 



For fc = 1 we have 



{q{i)i.v) - 



I 1^ 


-94 



93 





-94 


92 


-94 





-94 


aqi 


f793 





^^93 


92 



(2.25) 
(2.26) 
(2.27) 

(2.28) 



For higher k, the resulting lOM are the products of 9(1). 
For example, 

q(2)fj.ual3 = 9(l)aiy9(l)Ai/3- 

Thus, the higher Laxian lOM are not yielding any new 
independent integrals. 

An alternative variant of the Lax pair approach is pre- 
sented in Appendix H. 



D. I-Symmetry 

• Definition 1 ; The symmetry of the set of integrals 
{qa} (the I-symmetryJ is the mapping {x^} — > {x'^} 
which preserves the analytical structure {x'Y^^ = {{x'Y)* 
and leaves all the integrals invariant qa{x') — qa{x). 

I-symmetries define the algebraic structure of the sis- 
tem at hand. In practice one first calculates the infinites- 
imal I-symmetries, given by 



6qa{x) = Iqa = 6x^d^qa = 0, 



(2.29) 



where I = ijj^{x)d^ is the generating vector field, and then 
establishes the large (non-infinitesimal) I-symmetries, by 
exponentiating the infinitesimal ones. This is the stan- 
dard procedure in the theory of Lie-groups. 

Although in general the coefficients io^ are nonlinear 
functions of {x} (the nonlinear I-symmetries), here only 
the linear I-symmetry algebras will be considered. These 
can be calculated easily by a linear ansatz w'^ — a^^x'^ . 
In this way a general linear symmetry of the full set 
{qi, • • • , 54} of the regular lOM is found: 

25Bi = +^63^1 -I- (62 + iei)S3, 

25B2 = -^63^2 + (62 - iei)S4, 

2(5i?3 = -ie^B:^- a{e2- iei)Bi, 

26B4 = +ie3B4-aie2+i£i)B2. (2.30) 

In the spinor notation the matrix form of I-symmetries 

is 



4 



(2.31) 



where the traceless matrix S is given by 



1 



2 V £2 + »ei 



^(T(e2 - «ei) 



^leaSa, (2.32) 



and {ca} are real parameters. The basis matrices (the 
generators) are 



Si 

S2 

S3 



1 

2 
1 

2 
1 

2^3. 



cr 

1 

ni 

-i 



ifji forTG, 
-^0-2 forRG, 



ifTz for TG, 
for RG, 



2 



(2.33) 



The set of all matrices S that are traceless and satisfy 
the generalized hermiticity condition S^jy+TyS = forms 
the Lie-algebra su{2) for TG, and sm(1,1) for RG. The 
matrix 77 = ni is called the su(2)/ su{l, I) metric matrix. 
The generators {Sa} obey the standard commutation re- 
lations 



[81,82] = 
[82,83] = 
[83,81] - 



-iSi , 
-«82 . 



(2.34) 



Thus, both \ipi) and IV'2) are transforming according to 
the fundamental (spinorial) representation of the corre- 
sponding algebra gj. 

Every finite Lie-algebra g has the corresponding Lie- 
group G of "large" transformations, obtained via expo- 
nential mapping: 

\fSe g^G ■= exp(iS) e G. 

For gi = su(2) (the TG case), the group is Gi = SU{2), 
and for gj = su{l, I) (the RG case), the group is SU{1, 1), 
the noncompact version of SU{2). Both groups can 
be represented by sets of 2 x 2 complex matrices G 
that are unimodular (det G = 1) and (pscudo)unitary 
(GtryG^,?). 

The Cayley-Klein parameterization of the general 
SU{2)/SU{1, 1) group element 



G 



2/1 + m 
-o-(y3 - m) 



yi 



iy2 



(2.35) 



where yi,....4 G M, turns the unimodality condition into a 
geometric relation (the definition of the parameter man- 
ifold of the group Gi)- 



detG = {vif + [y2f +<j[{y:if + {v^f 



-- 1. 
(2.36) 



Thus, the parameter manifold for SU(2)i is the sphere 
S^, and for SU{1, 1)/ it is the hyperboloid H^, both em- 
bedded in M."^. (For a short classification of hyperboloids 
in see Appendix |c[) 



From this fact alone, one could expect that the TG case 
will be expressed in a natural way in terms of the trigono- 
metric functions, and the RG case in terms of both the 
trigonometric functions (compact dimensions) and the 
hyperbolic functions (noncompact dimensions). In this 
sense the cases arc "twins", i.e. there is a number of 
equations holding in both cases, up to the exchange of 
the trigonometric/hyperbolic functions. 



E. Action of I- Symmetries on the Lax variables 



It is of interest to know the action of the I-symmetries 
on the Lax matrices C^^i, . Since the Lax matrices are con- 
structed out of the basic spinors, some regularity must 
be induced in the transformation law of these variables. 

For example, for £11 = |V'i)(V'i|ni, the action of the 
infinitesimal I-symmetry yields 

SCn = S^-Cii +/:iinj;^S*ni. 

Owing to matrix identities rij^^ = rii and niS*ni = 
— S^, this expression simplifies to 



SCn = [S^.Ai] . 



(2.37) 



This represents the adjoint action of the I-symmetry on 
£11. In a similar way, one finds that the same transfor- 
mation law is valid for all C^^i,: 



(2.38) 



Thus, due to cyclic invariance of the matrix trace oper- 
ation, all Laxian lOM are invariant upon the action of 
I-symmetries. This is expected. 



F. Exceptional lOM 



The "regular" lOM, obtained in the subsection II B, 
form the full set of lOM for the complex coupling F. How- 
ever, in the special cases when F is either real or imagi- 
nary, there exist additional lOM. These will be called the 
"exceptional" lOM. 

To see the significance of these special cases, let us 
evaluate the "time" -change of the grating amplitude Q: 



IQ' = -TQq5. 

Here 55 is the expression 55 
"time" -change is 

Iqi, =A<jRe (F)|gp. 



(2.39) 

^1 + ^2 - (j(h + h), whose 



(2.40) 



Notice that is lOM in the case of imaginary F (so, it 
is an "exceptional" lOM). However, when F is a complex 
number, this quantity turns out to be a suitable variable 
for later calculations. 

.39|) and (2.4C) another conserved 



From equations 
quantity (for the general, complex F) is obtained: 



5 



q = ql + MQ\' 



(2.41) 



This quantity is lOM of the fourth order. It does not 
carry any independent information, since it is reducible 
to the already known regular lOM: 



(2.42) 



Nevertheless, plays an important role in one of the two 
presented procedures for solving EOM in the F S E case. 
From Eq. ( 2.39 ), two important relations follow: 



I\Q\' = -Re (F)|g|g5, 
/arg(Q)' = -Im(F)<75. 



(2.43) 



These equations indicate the existence of two important 
special cases: F S R and F G iM.. The case Im F = 
implies argQ = const, so that (/) := arg/i = const, while 
the case Re F = implies that |Q| = const. The F e K 
case is considered first. 



1. r e R 

Analysis of this case is based on the fact that the phase 
(j) of the grating amplitude Q is constant for real cou- 
plings. This allows introduction of a new independent 
variable 9{z) = dz'\fi{z')\ + 9q, which casts the prob- 
lem into a linear form. The matrix m is now replaced by 



m 







(2.44) 



const. 



where v := exp{iq 

The defining relation ( 2.11 ) has as solutions not only 
the matrix rii, but also the new one: 



I' 

-V 



(2.45) 



which is anti-Hermitiannl = —to-j,. Along the same lines, 



the defining relation (2.12) has as solutions both the ma- 
trix n2 and the new one: 



n4 



D 
av 



(2.46) 



which is symmetric. Having the new matrices na and 
n4 that satisfy the Lemma, a set of additional conserved 
quantities can be constructed: 



Wl 


EE (V'l 




= 2ilm {vBiBz) 


W2 


^ (V'2 


na \ip2) 


= -2ialm {vBiB2), 


W3 


= (V'l 


na \ip2) 


= —1JVB1B2 — 9B3B4 


W4 


= (^1 


04 iV'i) 


= vBl + avBl, 




= (^1 


n4 IV'2) 


= DBiBi-iyB3B2, 


We 


= {^2 


n4 IV'2) 


= vBl + avBl. 



(2.47) 



Note that {'ip2 \ na IV'i) = —w^, and (V'2| n4 IV'i) = ^5- 



One can extend the Lax procedure, in the spirit of the 



subsection [IC, to this case as well. The Lax matrices 
are now := l-^p )( V^i^l na, and the corresponding Lax 
equations 

9«/:« = [m, (2.48) 
The corresponding Laxian lOM of the first order are 



Tr£, 



-W5 



-W3 

W2 

-We 



W4 
W5 
—awi 
aw3 




(2.49) 



These lOM are the same as the ones already obtained 
through the spinorial approach. As mentioned, a more 
general Lax pair procedure is presented in Appendix 
Action of an I-symmetry on produces 



^ y;T/.(K) 



'-tip "3 



na. 



Here the condition n^'^H'^^n^ = — S"^, necessary for the 
covariant form of action, can be achieved in different 
ways: 

• Case 1: = —a. This corresponds to (j) = ±7r/2 (in 
TG) and to = or tt (in RG). This case allows for the 
full sui symmetry, i.e. all three Ca parameters can have 
non-zero values. However, only the diagonal part of Ti^ 
figures in the transformation law: 



(2.50) 



• Case 2: 7^ —a. Here only the diagonal part of 
Fsymmetries survives, i.e. ei and €2 have to be set equal 
to zero. The transformation law still has the same form 
as in the case above. 

Thus, the w lOM are invariant under the full sui sym- 
metry algebra if = —a, and under the ui subalgebra 
generated by cra/2 if v"^ ^ —a. 

An important special case is the phase conjugation, 
when the relative phase $ (= /3i -I- /32 — 71 — 72) is 
constant (0 or tt). Then, using relations (2.18) and the 
fact that the argument (= /3i — 71 = 72 — /32) of /i 
is constant, the following values for the integrals w are 
obtained: 

Wl = 0, 

W2 = 0, 

W3 = -y/qiq2c{a, ai - 02) exp(i(/32 - 7i))i 

W4 = qi exp(i(/3i + 71)), 

W5 = cr (j2 exp(?(/32 -I- 72)), 

we = -y/qiq2s{a,ai - a2)e-Kp{i{Pi + (32)). 



(2.51) 

These relations imply that all the phases /3i, P2, 7ii 72 
are constant, and that the a- variables are linearly depen- 
dent: «! — a2 = constant. Hence, all the fields essentially 
depend on only one real quantity, for example on ai. 
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2. r e iR 

The case of F imaginary has only one exceptional in- 
tegral, The corresponding I-symmetry of the set 
{qi , . . . , (75 } is restricted to the diagonal part of Eq. (2.30): 



2SBi 


= if^aBi, 


2SB2 


= -163-82 


2SB3 


= -iesBs 


25Bi 


= ^63^4 . 



This is a u{l) algebra of the transformations: 

S\lpl,2) = iyO-3 |'/'1,2) , 



(2.52) 



(2.53) 



and the corresponding group is U{1). The parameter 
space of this group is the circle §^ of circumference 47r. 
That group is the subgroup of both SU (2)/ and SU{1, 1)/ 
groups. 



III. SYMMETRIES OF THE EQUATIONS OF 
MOTION 



In general, one should distinguish the symmetries of 
the integrals of motion from the symmetries of the equa- 
tions of motion. 

• Definition 2 ; E-symmetries Any vector-field 
L that satisfies the master equation 



L.F 



0, 



(3.1) 



is the symmetry of the dynamical equations (I.4) 



The set gE of E-symmetries is also a Lie-algebra, i.e. it 
is linear, and the commutator of any two E-symmetries 
is another E-symmetry. So, one can describe the full al- 
gebra by its generators and their commutation relations. 

The E-symmetries are sought in the form of the most 
general linear ansatz 

L = x^a^^S^. (3.2) 

After some algebra, six generators are found for the 4WM 



system: 










io 


= -Bi5i + 


B2d2 + 


S353 + 


Bidi + c.c, 


Li 


= i{Bidi 


+ B2d2 


+ B3d3 


+ Bidi - c.c). 


L2 


- i{Bidi 


-B2d2 


- c.c), 






= ^(-8393 


-BaBa 


- c.c), 






= B2di- 


Bid2 + 


Bids- 


Bzdi + c.c. 


U 


= i{B2di 


-Bid2 


+ Bid3 


-Bsdi-c.c). (3.3) 



These six generators form the complete set of linear E- 
symmetries for the general (complex) F. They generate 



the algebra gs = r (B u{l) ® u(l) ® su(l, 1), with the 
non-vanishing commutators 

[-^1,-^4] = — 2-Z^5, 
[-^11-^5] = 2L4, 

[£4,-^5] = 2Li. (3.4) 
Defining the general infinitesimal E-symmetry by 5 = 
Ef=o ^i^i, we find the transformation law of the fields: 



S-B4 



[^0 
[^0 



'1 -r- 
h- 
h + 
h- 



]Bi 

]B2 
]B3 

IB4 



{O4 
[Oi 



or, in a more compact notation 



SIV'2) = 



% + 



iO^)B2, 
i9^)B^, 
i05)Bi, 

l05)B3, 



1^1 



(3.5) 



-*05)|^2 

l)1- j 







1^2) + 





-a{ei + i05)\A)- (3.6) 

The parameters and 0^ correspond to the two noncom- 
pact dimensions of the symmetry group Qe, i-e. their val- 
ues are arbitrary real numbers. This is in contrast to the 
rest of the parameters, which are periodic. So, the group 
of E-symmetries Qe = exp 5^ is M U{1)^ (g) SU{1, 1). 

It is easy to check [from the defining relation ( |3.lD ] that 
the E-symmetries always map the integrals of motion to 
the integrals of motion (and also the solutions of EOM 
to the solutions of EOM). Hence: 

Sgi = 2eoqi + (6*4 + ^6*5)94 + (6*4 - i05)qi, 
iq2 = 29oq2 - {Oa + i9^)qA - (6*4 - 16*5)94, 

S<73 = [200 + «(^?2 + 03)193, 

% - 2(0o + «0i)g4 + (04 + «05)(<Z2-gi). (3.7) 

This is a version of the Nother theorem: If one of lOM 
is taken as the " Hamiltonian" H of the system, then the 
action Li (H) of each i,; on such a Hamiltonian produces 
another lOM. 

Here, the following linear combinations of regular 
lOM are forming the irreducible representations of the 
u{1)lo'S)su{1, 1)li,L4,L5 algebra under the E-symmetries: 
qi + q2, qs and 93 form singlets 

^(gi + 92) = 20o{qi + (72), 

^qs = [200 + 1(92 " " 
5^3 = [20o-t{92 



while |T) 



^3)m, 

93)193, 



(3.8) 



{94, (91 ~ Q2) / V2, qi} is transforming as 
the triplet representation ^\T) — P\T), where P is given 
by 

-V2i9i + i95) 
29o 



( 2(00 + z0i) 
V2{9i-i95) 











-^/2{9i 



105 



V2{9i + i95) 
2(00 - i0i) 



Hence, one can start from the knowledge of only qi 
and recover the (almost) full set of regular integrals 
{<?i, 92, 94, 94}, by acting on them with the E-symmetries. 
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A. r G R 

In the r e R case, one can apply the same hnear ansatz 
as in the general case. The set of linear E-symmetries 
thus obtained is {Lq, • • • , L^}, where Lq-s a-re the already 
known generators (3.3), and the two additional genera- 
tors are found for the RG case: 

Lq = B^di- 3382 + 8283- Bidi + c.c, 
Lr = i{B4di +3382 + 8283 + Bid4-c.c.). 

(3.9) 

An alternative approach is to perform the redefinition 
z 6{z) of the "time" variable (introduced in the previ- 

av^ 
-V 

stant. This allows one to translate the vector- field lan- 
guage (applicable in the general case) into the matrix 
language. Define the column 



ous section), making the matrix rh 



I*) ^ I IS 

■\4'2 



(3.10) 



The evolution equation (2.7) can now be written as 



dg\-^) = M|*>, 
where the constant evolution matrix is 



M 




(3.11) 



(3.12) 



The master equation (3.1) now has the matrix form 

[K,M]=0, (3.13) 

where the matrix K = (K^j^) defines the infinitesi- 
mal symmetry of the big "spinor" i\^>) = Kj^*) (here 
H,v & {1,2,1,2}). The above master equation is trans- 
lated into "smaller" versions [K^^,, rh] = 0, valid for each 
2x2 block matrix K^^. Solutions of these "small" master 
equations are all of the same form 

K^t. = a^,^l2 + /3^i.m, (V/i, v). 



(3.14) 



The analiticity conditions 
constraints 



-era. 



a-i 



Pi, J ~ Pij I 



nin2 yield the 



(3.15) 



I.e. 



= -aafjl2 - o-Ajrii, 

Hence: 



(3.16) 



(3.17) 



In this way the rescaled EOM have 32 symmetries, char- 
acterized by the the real and the imaginary parts of the 
parameters {ofi^ , , a^j, = 1,2). 



B. On the relation between I-symmetries and 
E-symmetries 

One may ask the question, what is the relation be- 
tween the two groups of symmetries: I-symmetries and 
E-symmetries? The following general consideration clar- 
ifies this issue a bit. Let F be the EOM vector field, 
6 an arbitrary E-symmetry, § an arbitrary I-symmetry, 
and q an arbitrary lOM. Since [F,'S\ = and F{q) = 0, 
it follows that F{^q) = 0, i.e. fiq ^ q (Nother theo- 
rem: E-symmetry of lOM is also lOM). From this con- 
clusion and from 5q = Q it follows that [5, (5](7 = 0, i.e. 
[5, S\ ^ 5 (E-symmetry maps an I-symmetry into another 
I-symmetry). Hence, one expects that [Li,Sa] ^ St- 

This can be explicitly checked in the case of 4WM sys- 
tem: the generators Li for i € {0, 1,4,5} commute with 
all three Sa _generators, whereas the remaining two E- 
symmetries ^2,3 have nontrivial commutators with Sa- 



[L2,Si±S2] = T*(^l±^2j , 

[£3, 5^1 ±5^2] - ±1(5^1 ±52) , 



[1^2,3, S3] 



0. 



(3.18) 



Thus the su(l, 1) e symmetry (generated by {Li, L4, L^}) 
commutes with the su{2)j /su(l, l)j symmetry. 



IV. SOLUTION PROCEDURES 
A. r G R: The first procedure 

We present in detail the solution procedures for the 
case when F is real. This case is physically the most rele- 
vant. The case when F is imaginary, is treated similarly. 

Th e equations for ai and 02, extracted from Eqs. 
(1.4), form a closed system of equations: 



2Ia[ = 
2Ia'2 = 



-F [qis{a, 2ai) + q2s{a, 202) cos(<I>)] , 
-F [qis{a, 2ai) cos($) -I- q2s{a, 2a2)] , 



(4.1) 



which can be integrated with little difficulty. Once ai 
and a2 are known, the remaining four angles are found 
easily: 

2I(3[ = -arq2Sm{^)s{a,2a2)ticr,ai), 
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2/7; 
2/72 



(jTqi sin($)s(cr, 2ai)t{a, 02), 
-rq2 sin($)s(CT, 2a2)ct{a, ai), 
-Tqi sin($)s((T, 2ai)ct{a, 02), 



(4.2) 



where t and ci are the remaining two trigonomet- 
ric/hyperbohc functions, formed by using the rule (A4) 
(see Appendix A). 

Equations ( [4.l[ ) are integrated as foUows. First, two 
new variables are introduced: 



xi — c{a, 2ai), X2 — c(cr, 2a2). 

In terms of these variables — qiXi 
(4.1) become 



Ix[ = T[qi+ q2c{a, p) - Xiq^] 
Ix'2 = r [q2+ qic{a, p) - X2qB] 



(4.3) 

g2a;2 , and Eqs. 

(4.4) 



Note that, due to symmetry, only one of Eqs. (44) is 
independent. The solution of the other is obtained from 
the solution of the first one by interchanging qi and q2 ■ 
This, however, holds only w hen T is real. On the other 
hand, using Eq. ( ^.4lD , Eq. ( 2.4C ) can be written as: 



Ic^^^T{q-ql) 



(4.5) 



The integration of this equation depends on the geome- 
try. For TG, the total intensity is constant (/ = qi +92), 
so that 



dq5 



Tz 

T' 



For RG, the intensity is not constant, and 
qbdq^ 



q 



2 
95 



= Tz. 



(4.6) 



(4.7) 



The value of the integral in TG depends on whether q is 
larger or smaller than . For the first case: 



95 (-2) = Vq tanh tanh 
whereas for the second: 
95 (-2) = cotanh 



qm\ , 



Tz 



(4.8) 



cotanh 



qm \ , Vq 
V9 



Tz 



(4.9) 



In RG: 



ql{z) = ql{Q)eM-'^^z) + q[l-ex:p{-2Tz)]. 



(4.10) 



Once (75 is determined, Eqs. (4.4) for xi and X2 (i.e. ai 
and Q!2) can be integrated. The problem, therefore, can 
be reduced to the determination of one variable. Other 
variables can be solved in quadratures. To complete the 
solution, it remains to fit boundary conditions. This 
problem, however, is more conveniently addressed by an 
alternative solution procedure. 



B. r G K: The second procedure 

Another convenient method for solution of 4WM equa- 
tions is based on the linearization procedure (the replace- 
ment of the "time" variable z by the variable 9{z) = 



|/^(z')|dz'-|-6'o). Then (2.7) remains the same, but the 
matrix m — )■ m becomes constant [p, ^ v = cxp(i0)]. 



The explicit solution of Eqs. (2/7) is now 



c(ct, O) avs(a,Q) 
-vs{a,Q) c(tT, 9) 



(4.11) 



where Q = 9 — 9^^ and is to be determined from the 
boundary conditions. The subscript stands for the 
quantities evaluated at z = 0. The matrix in Eq. (4.11) 
explicitly displays the SU nature of the symmetry of so- 
lutions, and allows for an easy identification of Euler an- 
gles for the problem: a = (j>, f3 = 29, 7 = —(j>. In this 
formulation (real F) only one independent variable (9) is 
found necessary. The angle (j) is fixed by the boundary 
conditions. 

The evaluation of is facilitated by writing \Q\ and 
in terms of 9: 

\Q\ = ^s{a,29)/2, q^ ^ -<j^c{a,29). (4.12) 

The form of the solution is different in the two geome- 
tries, since / is constant in TG, whereas it is not in RG. 
At this point the symmetry in treating the two geome- 
tries is broken. The solution of Eq. (4.5) is 



tan(0) = tan(6'o) exp {^/qTz/I) , 
in TG, and 

sinh(26') = sinh(26'o) exp(rz), 



(4.13) 



(4.14) 



in RG. The procedure for evaluation of 6*0 is also different 
in the two geometries. We first present the TG case. 

The angle 9o is found when boundary conditions are 
applied to the solution given by Eq. ( 4.13 ). The condi- 
tions are that the fields are given at the opposite faces of 
the crystal: Aj{z = or z = d) ^ Cj. In OPC C3 = 0. 
Using these conditions, a number of auxiliary quantities 
is defined: 

U = |C4p-|Ci|2 + |C2p, 
V = |C4p_-|Cip-|C2p, 

p — 2C1C4 exp(— z0), 

a = exp(-V^d//), (4.15) 
(all real), and a shorthand notation is introduced: 

x^tSin{9d-9o), y^ta.n{9d + 9o). (4.16) 
There exists a rational relation connecting x and y: 

(4.17) 



ux + p vy — p 



V — px u + py 
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It is seen that x and y depend on and [through Eq. 
(4.13)] on q. However, there also exists a relation ex- 
pressing X (and likewise y) only through q: 



± 



q — a 



q — a 



q-b 
q — a 



1/2 



(4.18) 



where a = + , b = p"^ + v'^ , c = p{u — v) — 2p|C2p 
This relation is used to write an implicit equation for q: 



(4.19) 



where ^ = (1 — a)/2x, rj = (1 + a)/2y. Thus, given the 
boundary conditions, Eq. ( 4.19 ) is to be solved numeri- 
cally, to determine q. Given q, x and y are found, and 6q 
evaluated from the relation 



tanf 



(4.20) 



This completes the TG procedure. 

For RG, one finds two expressions for the modulus of 
the grating IQj at z = 0: 



IQol = tanh(erf)|qV(e-l) 



= Hsech(|Ci|2-MC4 



(4.21) 



where \C\^ = ^ |C, 



exp(rd), and now \p\ 



IC2C4I. This yields an expression for sinh(8(;): 

\p\ie-l) 



sinh(6d 



e(|CiP + |C4P) + |C2|^ 



(4.22) 



Using Eqs. (4.21) and (4.22), an expression for tanh(26'o) 
is obtained: 



tanh(26'o) 



sinh(2ed) 



e — cosh(20c;) 
This completes the RG procedure. 

c. r eiR 



(4.23) 



It i s useful to note that beside the equations ( 2.39 ) and 
( ^.43 ), one can derive an equation for the intensity 



//' = 2(f7 - l)Re {r)\Qf 



(4.24) 



From these equations it follows that in the F imaginary 
case a number of additional quantities is constant: \Q\, 
(75, /. The equation for the phase can be recast as 



/9^arg(0) = -Fgs, 
where F =: iF, and solved: 



arg(Q(z)) = arg(Q(0)) - t^z 



(4.25) 



(4.26) 



Since \^i{z)\ = \rQ{z)\/I{z) = |F||g|// is constant, one 
obtains an explicit expression 



H{z) = \fio\ exp(i0o - i^z), (4.27) 

where \^iol ■= \T\\Q\/I, 0o ^2 + argf + arg(Q(0)) 
and n := Tq^/I. 

Now the evolution matrix 1x1(2:) from the spinor EOM 
( ^.7| ) has the form 



m(z) 











and the formal solutions of EOM (2.7) are 

\Mz)) = uiz)\M0)). 



, (4.28) 



(4.29) 



The Ufz) matrix is the ordered exponential (see Ap- 
pendix y for a discussion) of m(z): 



U(z) 



exp dz'm.{z') 



(4.30) 



where the plus subscript indicates the path-ordered na- 
ture of the exponential. In practice, to obtain the explicit 
form of U(z) in terms of non-ordered quantities, one has 
to solve the initial value problem (IVP) 



a.u(z) 

U(0) 



m(z)U(z), 
1. 



(4.31) 



For the specific m(z) the explicit solution to this IVP 
is found (see Appendix ^) : 



U21 
U22 



expl-t— j ^cos(^- 

2aQ ( Q.z\ 
i — ^ exp — i — sm 

2crQ / Viz 



= exp 



.Viz 



cos 



T 

Sz 

T 
VI . 

— i — sm 



.^z 



4a|Q|2 



(4.32) 
It is easy to 



where S := Ty/q/I and q 
check that detU(z) = 1. 

In this manner, for known initial values \^i{0)), the full 
solution at later "times" z > is given by Eq. (4.29). 
However, by the nature of the 4WM system, one knows 
only the part of initial conditions. The system represents 
a split boundary value problem. 

• For the TG case the beams i?i = Ai and B3 = A4 
are entering the crystal sample from the z = side, while 
other two beams B2 = A2 and B4 = A3 are coming from 
the z = d side. Thus, only 1-01) is determined at z = 0, 
while \ijj2) has a fixed value at z = d boundary: 



IV'i(o))= ( ), \Md))^ 



-<jC2 



In this way: 

\Mz)) 



U(z)|Vi(0)), 
U(z)U(d)t|V'2(rf)). 



(4.33) 



(4.34) 
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• For the RG case only Bi = A\ and = ^4 are 
known at the z = boundary: -Bi(O) = Ci, i?4(0) = C4, 
while the remaining two field variables are given at the 
z — d boundary: B^id) = C2, Bj,{d) = C3. This means 
that both spinors j-i/'i) (« — 1, 2) are satisfying the mixed 
boundary conditions, where one component satisfies the 
initial value condition (at z = 0) and the other compo- 
nent satisfies the final boundary condition (at z = d): 



-aB2(0) 



\Md)) 
\Md)) 



Bi{d) 

C3 
B4{d) 
-aC2 



(4.35) 



where i?3(0), Bi{d), -82(0) and Bi{d) are unknown. In 
order to determine them, one has to use the evolution 



formula Eq. (4.29) to express the unknown boundary 
values in terms of the known ones. After some simple 
algebra one obtains 

Bi{d) 

^3(0) 





+ Ui2{d)C3 




U22{d) ' 


C3 


-U2l{d)Ci 




U22{d) 




-nUi2{d)C2 




U22(d) 


C2 


+ (jU2i{d)Ci 


U22{d) 



Here the unimodality condition detU(z) = 1 (Vz) was 
used. This concludes the solution procedure for both ge- 
ometries. 



symplectic structure on V^. In the phase space with the 
canonical coordinates {q^ , ...q^ ,pi, ...po} (9* not to be 
confused with the conserved quantities), the canonical 
form of F is 



F = ^ dgfc A dpk 



(5.1) 



fc=i 



The 2-form F sets up an isomorphism between the tan- 
gent space TV and the cotangent space *TV. Denote the 
inverse mapping by J : *TV — s- TV. In the canonical co- 
ordinates J has the form 



D 

fe=i 



(5.2) 



In a general system of coordinates {x^\fM = 1, 2D}, the 
forms of F and J become 



F = ^F^.dx'^ A dx". 



(5.3) 



with mutually inverse skew-symmetric matrices F = 
{Ffiiy) and J — [J^'^). Here the summation over the re- 
peated greek indices is assumed. The matrix J is known 
as the symplectic matrix. 

Physical quantities are smooth functions on V , form- 
ing the space C°°{V). A Poisson bracket is defined in 
C°°(y), generating a Lie-algebra structure 



{f,9}pB = J'"'d^fd,g. 
In the canonical coordinates this means 



(5.4) 



THE (PSEUDO)HAMILTONIAN 
STRUCTURE OF 4WM 



Considering again the form of the spinorial EOM ( p.7[ ), 
one may notice that (in the TG with /i e K. case) the ma- 
trix m is antisymmetric, resembling the symplectic ma- 
trix used in the Hamiltonian formalism for mechanical 
systems. This notice gives rise to the question whether 
it is posible to reformulate the 4WM system as a Hamil- 
tonian system. In this section one possible approach to 
the problem is considered. First the neccessary general 
definitions are given, and then the specifics of the 4WM 
system are discussed. 



A. Preliminaries 

For the formulation of the Hamiltonian formalism p!o| ] 
one needs a phase space in the form of a smooth mani- 
fold V and a closed nondegenerate differential 2-form F 
(the field- strength form) defined on it, which endows a 



{L9}PB = ^{dqkfdp^g - dgkgdpj). 



(5.5) 



fc=i 



The Poisson bracket is bilinear, skew-symmetric, and 
obeys the Jacobi identity, which is equivalent to the close- 
ness of the 2-form F: dF = 0. Later more will be elabo- 
rated on this condition. 

The dynamics is determined by the choice of the 
Hamilton function H on V. The external differential dH 
is a covector field (1-form), and J • dH is the correspond- 
ing Hamilton's vector field on V. The Hamilton equation 
of motion is specified by equating the tangent vector field 



X := x^dfj, with the Hamilton's vector field: 



= J -dH. 



(5.6) 



The Poisson bracket {/, g} pb may now be represented 
by the action of the covector df on Hamilton's vector 
field J • dg: {/, g}pB = df{J ■ dg). Therefore, the deriva- 
tive of function / in the direction of Hamilton's vector 
field J ■ dH is in fact {F,H}pb- Hence, the Hamil- 
ton equation ( ^.6| ) can be written as / = {f,H}pB for 
an arbitrary function /. Since the coordinate functions 
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{q^, ...q^ ,Pi, ...pd} form a complete basis, the equations 



g^- = {q^H}pB^dp,H, 
Pk = {pk,H}pB = -d„kH, 



(5.7) 



form a closed system. These are the canonical Hamilton 
equations of motion. 



B. Four- wave mixing 

In an attempt to cast the 4WM system in the Hamil- 
tonian form, one encounters several problems. 

First, there is no clear choice of the Hamiltonian 
(Hamilton's function) H{x). It can be an arbitrary real 
function of the full set of conserved quantities: H{x) = 
h{q{x)). For example, for the general F G C, one can 
identify three convenient families of Hamiltonians 



Hqi^e = qi+ £92, 
Hq^\ = Xqs + Xqs , 

Hq4X = Mi + Mil 



(5.8) 



where e = ±1 and |A| — 1. Clearly, these families are 
not exhausting all the possible choices, even among the 
Hamiltonians that are linear in the regular lOM. 
The corresponding symplectic matrices are 



'g3A — 



(5.9) 



-Xfi 



Xjj, 



aXfi 



>q4,\ 



-Xjj, 



A/i 



—A/i 



-A/t 



A/t 



A/i 



V 



-A/t 



-A/i 



The matrix elements that are not written explicitly, are 
zero. Note that for these three families (subscript K is 
the index of each individual family): 



'K 



_tt 



'K 



IK- 



(5.10) 



These matrices are chosen in such a way to satisfy the 
three basic requirements: they are antisymmetric, non- 
singular (in the matrix sense), and they give the same 
EOM 



J^'O^Hk = Fix), 

where f^ix) is the right-hand side of the EOM (0): 
x'^ = ff^. The fact that all Hamiltonian structures must 
reproduce the same dynamical equations ( [l.4[ ) means 
that for any two structures (Ha, 3 a) and (Hb,Jb) one 
can write: 



^WbaHa = VHe 



(5.11) 



where '''Vba ■= ^ba • V defines the cogradient and 
TIba '■= J^^Ja is the recursion matrix, connecting the 
Hamiltonian struc tures {A) and {B). 

Equation ( 5.11 ) defines the mapping from the first 
Hamiltonian structure {A) to the second one {B). It 
is interesting to assume for a moment that there exists 
some function He whose gradient is the BA-cogradient 
of Hb- If such a function exists, it will be an integral of 
motion: 



d.Hr 



x^'d^Hc - 

j^-'dMJB': 

a/3; 



f^aJA^dpHB 



— —daHB-JA dfjHB — 



Thus, the conserved quantity He, if it exists, is a new 
Hamiltonian of the system, and the corresponding sym- 
plectic matrix is Jc ■= RabJb = JbJ^^Jb- 

One can continue along the same lines, defining a series 
of conserved quantities ("Hamiltonians") and the corre- 
sponding symplectic matrices: 



Ha 
Ja 



Hb 
Jb 



He 
Jc 



(5.12) 



The sequence terminates when the Hamiltonians start re- 
peating themselves (i.e. they became linear combinations 
of the previous members of the sequence). This type of 
sequence of the Hamiltonians and the corresponding sym- 
plectic matrices is common in the two-dimensional inte- 
grable systems. There exists a multi- Hamiltonian prop- 
erty of integrable systems, whereby the chain of Hamilto- 
nians is (usually) non-terminating, leading to an infinite 
set of non-equivalent lOM, and, thus, to the complete 
integrability of the system. 

In 4WM one expects that all such sequences, if they 
exist at all, should terminate after a few terms. For ex- 
ample, consider the two structures {Ha ■= -ffgi.+i = 
91+92, Ja) and {Hb := i?Qi,-i =^ gi -g2, Jb), where Ja 



and Jb are the special cases of (5.9). The recursive ma- 
trix is Kba ■= Jb^Ja = diag(l,-l, 1,-1, 1,-1, 1,-1) 
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and the basic identity is ^VbaHa = ^Hb- The defini- 
tion of the induced He is 



'baHb 



(5.13) 



and its solution is He — qi + q2- Thus He = Ha and 
the sequence is periodic: Ha Hb Ha 

Each Hamiltonian structure {Hk,Jk) has the corre- 
sponding Poisson bracket: 



{f,g}K :^ (WffjKi^g). 



(5.14) 



This bracket is antisymmetric {g,f}K = ~{f,g}K, and 
can be characterized by the set of basic brackets: 



K 



(5.15) 



However, a bracket so defined does not satisfy the ex- 
pected Jacobi identity {/, {g, h}K}K + cyclic = 0. The 
tensor of non-Jacobianity measures the violation of the 
Jacobi identity: 



eil^iva 



{{x^, x^jK, x^jK + cyclic 
{d^J^") J"" + cyclic. 



(5.16) 



This tensor is essentially the same as the tensor of de- 
flection from the Bianchi identity, to := dF (i.e. w^j/a ~ 
cyclic): 



8 F 



"(if) 



(5.17) 



In the case of non-singular J, the Jacobi condition is 
equivalent to the Bianchi identity. 

In the system at hand, the fact that ri(x) is not dis- 
appearing is the consequence of the non-constancy of fi. 
For example, take again the Hamiltonian H = i/gi.+i = 
qi + q2- Then 



(5.18) 



where 



1 \ 





-1 

V 



E_ 



\ 

1 



Note that E^ = -E±, E± • E^ = 0, and 

9 1 

(E±) =--[l8±fT3®CT3 0fT3]. 



(5.19) 



As a simple consequence of these expressions, the follow- 
ing identity is valid 



(aE+ + /3E_)"^ = --E+ - ^E 



(5.20) 



for arbitrary (nonzero) a and (3. This identity is used to 
evaluate the "field strength" matrix: 



F = --E+ - 4e_ 
From this, one obtains 



(5.21) 



= ^ [{dafi){E+)^^ + cyclic] - 



(5.22) 



For example, i^Bi, 34,82 ~ ~Jfi'^^l^- Thus, the Bianchi 
identity is clearly broken, and one can not find the po- 
tential Ap,{x) such that Ff^i,{x) is its strength tensor. 

One can see the dual nature of the same obstacle, ex- 
pressed in terms of the Poisson bracket, in the following 
way. The basic brackets are 



IJ', 
fl. 



and one non- vanishing component of O is 



(5.23) 



(5.24) 



The full list of non-vanishing components of is given 
in the Appendix ^. Thus, the Poisson bracket is not 
self-consistent: one can not apply it consecutively on the 
phase space without running into inconsistencies. This 
is the second, and much more serious problem with the 
presented approach to casting the 4WM system in the 
Hamiltonian form. One can say that the 4WM system 
has a pseudo-Hamiltonian structure. 

If all components of u were zero, one would be able 
to find the potential functions Af^ of strength tensor 
F^u = d^Au — di,A^. Then the system could be formu- 
lated as the Lagrangian one, with the action functional 
S{z\,Z2) := J^^ dzL{x,x), where the Lagrangian func- 
tion is 



L(x,i) := xf'A^ix) - H{x). 



(5.25) 



The Euler-Lagrange EOM corresponding to this La- 
grangian are the Hamilton equations ( |5.6| ). The elements 
of the Lagrangian formalism are provided in Appendix 

Since uj ^ 0, one may search for the solutions in the 
form 



F 



fid^A,~d,Af,), 



(5.26) 
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which leads to 

{daf)Ff_,„ + cyclic = fuap-i- (5-27) 
The direct consequence of this equation is 



da In/ 



1 



2D -2 



(5.28) 



where 2D = 8 in the 4WM system. After some algebra 
one derives 



i9a In / = -da In + ^ (^3 ® c^a ® (^?.)a "(^i^^ 



(5.29) 



where (j) = £^rg(/^)- 



C. The r G R case 



If (f) is constant (i.e. F e M), the solution of Eq. ( 5.29| ) 
is / = Then 



Note that if in the above expression Ai (z) is directly re- 
placed by |/;i(a;)|, the obtained corresponding variation 
equations are wrong. 

In the r G K case the set of lOM is enlarged by the " ex- 
ceptional" lOM {wi-e} (and their complex conjugates), 
and one can construct some additional families of (linear 
in lOM) Hamiltonians (e = ±1, |A| = 1, \e\ = 1): 

Hwi,€ — wi + ea_W2, 
Hw3,\ = + XW3, 

Hw5,X = AW5 + XW5, 

Hw4,xe = -{Xwi + Xwi + 6wq + 9wq) . (5.33) 
The corresponding symplectic matrices are 
/ 



where F^^ := F^^/ f — d^Au ~ d^^A^, and the solution 
for the potential A^ is 



An 



1 r 



(5.30) 



with v = exjp{i(j)) (not to be confused with the index v). 

To construct the action for this case one has to go 
one step back. The factorization of / = l/|/i(x)| from 
F^^ is equivalent to the introduction of a new time pa- 
rameter 9{z) := dz'A4{z') + const, into EOM, where 
M{z) ■= \^{x{z)) I is the on-shell value of |Ai(a;)|. Thus 

F,,^=d,H{x). 



Act 



Act 



Act 



-aX 



-A(TN 



The constant "field-strength" form F is closed, and its 
tensor of non-Jacobianity uo disappears. So, one can con- 
struct the action in the rescaled time 



^(^1,^2 



(5.31) 



Further transformation of this action 

fZ2 



s 



dzM{z 
1 



1 ^ ~ 1 dx" 
2^ ^"•'Mi^^'^^^'l 



dz 



dr^ 

x^F,„,^-M{z)Hix) 



dz 



leads to the Lagrangian 
1 



-M{z){qi+q2), 



(5.32) 



-A 



V 



J 



These matrices 3w--- do not depend on /i, i.e. they are 
constant. So, their Poisson brackets satisfy the Jacobi 
identity. This is the case not only for the real, but also 
for the complex coupling F. However, the Hamiltonians 
( p. 33 ) are not lOM for the complex F. 

To sum up the results, for a general F two types of 
pseudo-Haniitonian structures exist: 

• {Hq,3q): Hamiltonians Hq are linear in regular 
lOM. They are conserved quantities in general case, but 
the corresponding field-strength forms Fq are non-closed 
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dpQ ^ 0. The defect of this type is the non-closeness of 
its symplectic structure. 

• {Hw,-fw)- Hamiltonians Hw are hnear in excep- 
tional € K.) lOM. The corresponding field-strength 
forms Fw are constant and closed in general case. The 
defect of this type of structure is the nonconservation of 
VF-Hamiltonians (in general, F €E C case). 

In the case of real coupling, both defects disappear, 
the first one after rescaling z 0{z), and the second 
one because T^-Hamiltonians become constant. Then 
one can construct a consistent Hamiltonian structure for 
the 4WM system. 



VI. CONCLUSIONS 

In summary, the algebraic structures of the 4WM equa- 
tions in PR crystals were studied. 

First, the form of the equations of motion was used to 
group the basic fields into two doublets, leading to the 
new form of EOM, resembling the Dirac equation in one 
dimension ("time"), with the field-dependent mass ma- 
trix. This lead to the simple procedure for finding the 
complete set of "regular" (i.e. present in the complex F 
case) integrals of motion. Then an alternative but closely 
related procedure, based on the Lax pair approach, was 
used to check the completeness of the obtained set of 
lOM. Both procedures were extended to the special case 
F e M, to obtain an additional ("exceptional") set of 
lOM. 

Afterwards, the concept of symmetries of the "regular" 
lOM was defined (the I-symmetries) , and the Lie algebras 
(and groups) corresponding to the linear I-symmetries 
were found. These are the su{2) symmetry for the trans- 
mission gratings and the su{l, 1) symmetry for the reflec- 
tion gratings. The initial doublets of basic fields, which 
were introduced as a convenience for more compact cal- 
culations, turned out to be the fundamental (i.e. spinor) 
representation of those symmetry algebras. Also, the Lax 
matrices, constructed from these basic spinors, transform 
in the regular way, i.e. they form the adjoint representa- 
tion of the I-symmetries. 

In the special case F e M the number of lOM increases, 
so only the subset of "regular" I-symmetries survives. 
This is to be expected, since the I-symmetries now have 
to satisfy a larger set of constraints than in the general 
(F e C) case. 

In the second part of the paper another type of sym- 
metries was considered, the symmetries of EOM (the E- 
symmetries). The corresponding symmetry algebras are 
the products of several abelian factors (one noncompact 
^ M^, and two compact ~ u{l)) and of one su(l, 1) factor 
(for both geometries). The action of these symmetries on 
the regular lOM was studied and a special kind of Nother 
theorem is found to be valid here. 



In the special case F G R the number of independent 
EOM gets smaller, leading to the increase in the number 
of E-symmetries. This is clearly the opposite behavior 
to the case of I-symmetries. Further study is necessary 
to clarify the relation between the "regular" and "excep- 
tional" cases. At the end of this part, the action of the 
E-symmetries on the I-symmetries was considered (in the 
"regular" case). The non-abelian factor of E-symmetries 
commutes with the I-symmetries, and the two uil) fac- 
tors act as rotations in the 1 — 2 plane of I-symmetries. 

As a short excursion from the algebraic orientation of 
the paper. Section IV is devoted to the solutions of EOM 
in two "exceptional" cases: F g R and F e iM. In both 
cases it is relatively straightforward to obtain the gen- 
eral solutions (two methods for F G M were presented 
and one for F € iM), but satisfying the boundary con- 
ditions characteristic of 4WM geometries required more 
attention. 

In the last Section one possible approach to the Hamil- 
tonian formulation of the 4WM system was discussed. 
The problems that occurred in that program were two- 
fold: the non-uniqueness of the choice of the Hamiltonian 
(Hamilton's function) , and the non-closeness of the field- 
strength 2-form. The first problem leads to the recog- 
nition of the multi-Hamiltonian nature of the 4WM sys- 
tem, and is not really a problem. It is just the type of the 
"gauge-symmetry" of EOM. The second problem, how- 
ever, is the real obstacle to the fulfillment of the program. 
The structure of this obstacle is topological (the violation 
of the Jacobi and Bianchi identities). This was studied 
for one specific " gauge" (the choice of the Hamiltonian) , 
and a special circumstance when this obstacle can be re- 
moved was found, essentially corresponding to the F G M 
case. 

The same F e M case was then treated in a different 
way, leading to the discovery of even bigger space of pos- 
sible Hamiltonians. The topological obstacle is absent in 
this case, and nothing prevents a full consistent applica- 
tion of the Hamilton formalism, and identification of the 
corresponding Hamiltonian action of the system. 

Future work: The presented work contains several 
topics that deserve future attention. 

Questions pertaining to the general class of dynamic 
systems: Clarifying the freedom of choice of Hamilto- 
nian function among lOM; Studying properties of E- 
symmetries upon the local (i.e. x-dependent) scalings 
of the dynamic vector field F; Finding classes of equiva- 
lency of the symplectic form J (under such scalings) that 
have the same structure of the tensor of non-Jacobianity 
LO] etc. 

Questions related to the 4WM system in particular: 
Full relation between the F S M and F e C E-symmetries; 
Explicit resolution of boundary conditions in F e C case; 
Extending the theory to multiple gratings; etc. 

The last question is particularly intriguing. Even in 
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the case of single gratings, 4WM EOM possess rich al- 
gebraic structure. However, the writing of gratings in a 
photorefractive crystal is a dynamic holographic process, 
and more than one grating can coexist simultaneously in 
the same region of the crystal. EOM then contain terms 
coming from different types of gratings, and the analysis 
should be much more involved. 

Acknowledgements: One of the authors (PLS) ex- 
presses gratitude to the Brown University for support 
during the graduate years, when some of the ideas ex- 
plored in this work were conceived and partially devel- 
oped. 



with arbitrary complex numbers a, /3 and 7. Its evolution 

equation is 



d,\^)=Af\^), 



where 











(B2) 



(B3) 



The problem is to determine the evolving member of 
the Lax pair. 



APPENDIX A: THE a-METRIC ELEMENTARY 
FUNCTIONS. 

The elementary definitions and relations of the c and 
s functions are listed in this Appendix: 

cos(ga), 
cosh(Q:), 

sin(a), 
sinh(a), 

(Al) 



c(cr, x) := cos{^/ax) = 
s{a,x) := sm{-\/ax) / y/a 



where the upper/lower option corespond to a = +1/ — 1 
signs. 



c(cr, x)^ + as{a, x)"^ = 1, 

2c((7, x)s((T, x) = s((T, 2a;), 

2c{a,xf = l + c(a,2x), 

2as{a,x)'^ = 1-c(ct, 2a;), 

c(a-, a;)^ — C7s(cr, a;)^ = c(cr, 2a;), 

t((7, x) := s{a,x) / c{a,x), 
ct{a,x) := l/t{a,x), 

t((T, x) + (Tct((T, x) = 2a/s{a,2x), 
t{a, x) — act{a, x) = —2act{a,2x), 

c{a,x) = —as{a,x), 
s'{a,x) = c{a,x). 



(A2) 



(A3) 

(A4) 
(A5) 
(A6) 



APPENDIX B: AN ALTERNATIVE LAXIAN 
APPROACH 

In order to achieve a sufficient degree of generality, one 
should use the "big spinor" 



(Bl) 



1. £~ |*)(*| 

The matrix £ is searched first in the form £ = 
^ (\E'| H, where A and B are some constant matrices. 
Then 



Require 

B-^Af'^B = -AJ\fA-\ 



I.e. 



Af'C = -CAT, 

where C := BA. 

The solution of this equations is 



(B4) 

(B5) 
(B6) 



C = 



/ "6 
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(B7) 



In the general case of complex F the factor is non- 
constant (with respect to z) , and some of the parameters 
^ have to be set to zero: = ^3 = ^5 = ^7 = 0. However, 
if r is real, the phase factor v.= cxp(iarg/i) is constant, 
and one may redefine the time variable z ^ 6*, to absorb 
the non-constant absolute value |m('2)|- In effect this per- 
mits the full set of non-zero ^ in the above matrix C (with 
the replacement jj ^ v). 

Let us consider the general case. The presence of 
four non-zero ^ parameters indicates the existence of four 
lOM: 



JOMj, = a{Ii+ als) = aqi, 

lOM^^ = /3j{l2 + ah) = f3^q2, 

lOM^^ = aPiBiBi - B3B2) = ddq?,. 

lOM^^ = a^{BiBi - B2B3) = a^qs, 



where 



lOM,, :. Tr- ^ ( * 



dC 



(B8) 



(B9) 
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These integrals are the already known "regular" TOM. 
Only 54 is not obtained in this way. It will be obtained 
in the next subsection, with a different choice of L. 

In the r e M case, there are eight free parameters (^) 
and there should be eight lOM. The first four of them 
are the "regular" ones {(71,92,93,93}, and the additional 
four are: 

lOM^, = a\-r\H}^BiB2-PB4B2)^-h\^W2, 
lOM^^ = aj{vBiB2 + <TvB^B4) = —aajw^, 
lOM^^ = —aajw^. 

(BIO) 

Here, again, all obtained integrals are already known. 
They are the elements of the w-set derived in the spino- 
rial approach. The remaining elements of that set will be 
obtained in the next subsection. 



APPENDIX C: HYPERBOLOIDS IN 



In Mf^ there exist four different types of the normalized 
hypcrboloids" , defined by yf + £22/2 + ^32/1 + ^^vl = 1- 





£2 


£3 


£4 




+ 1 


+ 1 


+ 1 


^^3,1) 


+ 1 


+ 1 


-1 


^(2,2) 


+ 1 


-1 


-1 


^(1,3) 


-1 


-1 


-1 



^(4 0) j^^^ another name for the sphere and M'^^ 
is the hyperboloid relevant for this work. 



2. £ ~ I*) (*| 



APPENDIX D: ORDERED EXPONENTIAL 



Now search for £ in the form £ = ^ |^) (^| ;B, where 
A and B are some constant matrices (different from the 
ones in the previous subsection). Then 

d,£^AAfA-^C + CB-^Af^B, (Bll) 

and the requirement that C satisfies the Lax-type evolu- 
tion equation has the form: 

B-'^M'^B = -AJ\fA-~\ (B12) 

i.e. 

Af^C = -CM, (B13) 
where, again, C := BA. The solution to this condition is 



C = 



/ap.^2 


6 




ai5 \ 






(3 ?5 




13 ^6 










-6 





(B14) 



In the r e C case one has to set ^2 = = = ^4 = 
and the remaining ^ parameters give rise to the following 
four "regular" lOM: 



JOM53 

lOM^, 

The r e R 
lOM^^ 
lOM^, 
lOM^, 
lOM^, 



= 0, 
= 0, 

= aj{BiB2 + aB^iBi) = ajqi, 
= ajqi. 

case has additional w integrals: 
= a[DB\ + auB^) = awi, 
= f3"fia9Bl + vBl)= Pjwe, 

vB-iB2) = a(3w5, 



(B15) 



aPiPBiBi 

IW5. (B16) 

In this way, the full sets of "regular" (9-set) and "excep- 
tional" (ui-set) lOMs are reconstructed. The presence of 
the arbitrary complex constants a, /? and 7 in the pro- 
cedure indicates that there are no additional lOM of the 
bilinear type. 



In this appendix some general properties of the matrix 
U(z) are discussed. 



The basic 4WM EOM (|/7D has a formal solution 
\i,,{z)) - U(z) 1^,(0)) , (Dl) 
where U(z) satisfies the initial value problem 



9,U(z) = m(z)U(z), 
U(0) = 1. 



(D2) 



One can write the formal solution to this equation in the 
form 



U(z) := ( exp ( I dz'm(z') 



(D3) 



which is called the (Path) Ordered Exponential (OE). The 
notion of ordering is referring here to the right-to-left 
multiplication of the factors in the definition of OE: 



u 

j^g(J„^dz'm(z'))^ ■= ^im [| e(^'"("^)), (D4) 

a=N 

i.e. one alternates the infinitesimal integrations (along 
the path between the z' = and z' = z) and exponenti- 
ations of, in such a way obtained, infinitesimal matrices. 

OE is an entirely different object from the ordinary ma- 
trix exponential exp(Jp c?z'm(z')), where the whole inte- 
gration along the path is performed first, and then only 
one exponentiation executed on this integrated matrix. 
The source of the difference is in the non-commutativity 
of the matrices m(z) evaluated at different points. 
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stant matrix m(z) = 



and all commutators 



The methods to evaluate OE are frequently non-exact: 
one may easily prove that the knowledge of U(z) for arbi- 
trary m{z) is equivalent to the knowledge of the solution 
to the Schrodinger equation for an arbitrary complex po- 
tential (and this is known to be a non-solvable problem) . 
However, in the cases when m(z) has one of the several 
special forms, the exact solution for U(z) can be found. 
Two such cases are encountered in this work: 

• for r G M the matrix m(z) is proportional to the con- 
av" 
-V 

[m(z), m(z')] are equal to zero. Thus, OE reduces to the 
ordinary exponential, and the result is displayed in Eq. 

• for r e iM the matrix m(z) has the raising and the 
lowering components that oscillate with the opposite fre- 
quencies Q.. The Appendix g gives one possible way to 
obtain U(z) for such an m(z). 

In the general case, the OE U(z) satisfies several sim- 
ple identities, induced by the properties of m(z): 

• From the tracelessness of m.{z) follows the unimodal- 
ity condition detU(z) = 1 (Vz): 



detU = lim TT dct 



lim 



a=N 


n 

a=N 




exp 



= lim TT exp(O) = 1. 



(D5) 



• From the membership of m(z) in the Lie alge- 
bra g := su ((3 + (t)/2, (1 - o-)/2) it follows that U(z) 
is an element of the corresponding Lie group Q := 
SU{i3 + a)/2,{l-a)/2): 



(D6) 



where the 77 = ni^^^ ^/^^ metric matrix of the 

algebra g. To check this identity, one should follow the 
chain of arguments: 



m 



^ exp(^m(a^))eg 
^ n exp (^m (4)) eg 

Viz) eg. (D7) 



APPENDIX E: ORDERED EXPONENTIAL 
SOLUTION FOR T e IR 



In this appendix the initial value problem 

d,V{z) = m(z)U(z), 
U(0) = 1, 



(El) 



is solved for F imaginary. Starting with the ansatz 

Uij{z) = Mil exp{iujijz), (E2) 

one obtains a set of conditions 

CJ21 = iJ^ii + ^, 
UJ12 = LO22 — i^, 



U21 

Ull 

U12 



-Mo 



IUJ21 

afio 



(E3) 

U22 ~f^0 1^12 

The second pair of these conditions defines the consis- 
tency conditions 

2 



^^11^21 = LO12UJ22 = O-IjUo 

which are converted into auxiliary equations 



^11 



Slwii 

f2w22 



CTlMol 



= 0, 
= 0. 



(E4) 



(E5) 



U22 — iiu/22 — CT|^o| 

Solutions to these quadratic equations are 

wii± = cji2± = (-rj±s)/2, 

c^2i±=W22± - (+r!±S)/2, (E6) 
where S := T^/I and q := ql + 4cr|(3p. Then 

Utj = Uy+ e-xp{izijJij+) + Uij- exp(iza;.y_), (E7) 
and one has to solve the remaining conditions 

U21± iCTUJll± 



Ull± 
Ul2± 



Mo 
-«a;22± 



U22± Mo 
in conjunction with the initial conditions 

1, 
0, 
0, 



(E8) 



U11+ 
U12+ 
■"21+ 
■"22+ 



Ull- 
"12- 
U2I- 
"22- 



1. 



(E9) 



The solution 


is 






U11+ = - 
U22+ = - 


wii- 


"11- = 
"22- = 


^^11+ 


S ' 

ti^22- 


E 

^^22-1- 




E 


U12+ = - 


tCTflQ 


"12- = 




E ' 


E 


U21+ = 


«Mo 
^ E ' 


"21- = 


«Mo 
E ' 



(ElO) 



leading to the final form 
C/11 

U12 

U21 

U22 



exp — 

.2aQ 
i exp 

.2aQ 
-i — — exp 



cos 

'riz 
Y' 
riz 



n . 

i — sm 



exp 



nz 



cos 



^z 

T 



— i — sin 



(El] 



It is easy to check the unimodality condition det U(z) 
(Vz). 



18 



APPENDIX F: COMPONENTS OF Q 



The non-vanishing components of 17 are 





= — /i94/i, 




= — 




= +Aic^2A*) 




= —fldzn, 




— ~iidi^, 




= +fj,di^, 


QB1B3B4 


= +p.djifi, 


QB^BiBi 




^^2^4-63 


= +p.d_ip., 


^^2^4-64 


= +^J,d2fi, 


QB2-B4-B1 




^^2^4-62 


= -iidAp., 


QS2S4-B3 


= +fJ.difl, 


^-62 ^4 ^4 


= +p,d2p.. 



(Fl) 



APPENDIX G: LAGRANGE FORMULATION OF 
GENERALIZED HAMILTONIAN SYSTEMS 

1. From the singular Lagrangian to the generahzed 
Hamiltonian dynamics 



The singular Lagrangian (/i e 1,-0): 

L = A^ix)i''-Vix), (Gl) 
generates the Euler-Lagrange equations of motion: 

where := 9^.4,^ — 9^^^^ is the tensor of the "field 
strength" . If it is nonsingular (det F 7^ 0; possible only 
for even D), its inverse, the tensor of symplectic struc- 
ture J := F^^ can be defined. Then Eq. ( p2[ ) have the 
form of the generalized Hamilton equations: 



(G3) 



To initiate the Hamilton formulation, one starts with the 
definition of canonical momenta 



dL 

dxf^ 



and of the (naively defined) Hamiltonian function 
H := p^x^ - L. 



(G4) 



(G5) 



The definitions of momenta (G4) are supposed to be 
(non-singular) contact transformations replacing veloc- 
ities x^ by the momenta p,^, thus allowing (by inversion) 
to express the velocities x^ in terms of the momenta 
Pfj_. Here, however, these equations are singular: the ve- 
locities do not figure (at all) on their right-hand-sides 
(RHS). Thus, instead of being the (successful) contact 



transformations, these equations are the constraints on 
the Hamiltonian dynamics of system: 



:J3m -A^,{x). 



(G6) 



Constraints obtained from the contact transformations 
are called the primary constraints, implying that 
there may be some additional constraints in the system 
(all these additional constraints are called secondary 
constraints) . Hamiltonian systems with constraints are 
treated by the Dirac method which outlines are given in 
the rest of this Appendix. 

For the specific system at hand, the naively defined 



Hamiltonian (G5) has the arbitrary weighted terms pro- 
portional to primary constraints: 



H 



V 



(G7) 



Its unique (non-arbitrary) part is called the Canonical 
Hamiltonian He = V . For the further purposes, one 
needs the (temptative) Total Hamiltonian: 



H'rp — He 



(G8) 



where are the Lagrange multipliers, for this mo- 
ment taken to be arbitrary functions of time z. The 
prime on H!j, denotes the temptative nature of this quan- 
tity: once when (and if) the Lagrange multipliers are 
determined (i.e. replaced with suitable functions over 
phase space), this quantity will be replaced by symbol 
Ht- One can not know apriory what values will vs ac- 
quire. Instead, the (Dirac's) constraint analysis has to be 
performed to determine the full set of constraints in the 
system and then to classify the constraints as either the 
first class (ones that commute in a weak sense with all 
other constraints) or the second class (ones that are 
not of the first class). Then, the Lagrange multipliers 
standing next to the primary constraints of the second 
class will be determined as a specific functions on the 
phase space, while the Lagrange multipliers correspond- 
ing to the primary constraints of the first class will stay 
undetermined. Presence of the second class constraints 
leads to the reduction of the phase space, while the first 
class constraint generate the gauge transformations (on 
the phase space). 

The phase space F — {{x,p)\Vx,p G M^} possesses the 
Poisson bracket; 



{f,g}pB ■■= 



df dg df dg 



^x^^ dpa dpn dxt^ 



(G9) 



The evolution of any function / on the phase space F is 
defined by 



f^={f,H!,}pB. 



(GIO) 



The consistency equations of the primary constraints 
dKn 



dt 



^ = {k^, H^}pb = -d^V + F^.v'' = 0, (GU) 
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can be solved if F = is non-singular, giving the 

final expression for the Lagrange multipliers: 



(G12) 



I f F is singular, some components of the RHS of Eq. 
( |G11 ) can not be solved for v, and one has to define 
the secondary constrains (and then to check their con- 
sistency, and so on). In this work, only the case of non- 
singular F is discused. 

Here, all primary constrains are of the second class, i.e. 
each of them has (at least) one other constraint that does 
not commute (in Poisson bracket sense) with it. This is 
easy to see from 



(G13) 



and the non-singularity of the matrix F. Constraints 
of the second class reduce the phase space of the system 
(whereas the first class constraints, if they existed, would 
be the gauge symmetries of the corresponding action). 

To successfully perform the reduction of phase space, 
one needs to replace the Poisson brackets with the Dirac 
brackets, defined as: 



{f:9} 



DB 



PB 



(G14) 



2. When the generalized Hamiltonian dynamics has 
the Lagrangian formulation? 

One can turn any even-dimensional generalized Hamil- 
tonian system (given by Eqs ( |G3| )) into the Lagrangian 
(^ iff:' 

a) Matrix J is nonsingular, so one can define its inverse 
F, and 

b) 2-form F := ^F^i.dx'^ A dx" is closed, i.e. satisfies 
the Darboux condition dF = 0. 

The two conditions are (in a simple connected region 
U of the phase space: Tri{U) = 0) sufficient to assure the 
existence of "potentials" A = Afj.dx'^, such that F is ex- 
act 2-form F = dA (and, therefore, closed). Under these 
conditions the potentials are: 



A^,{x) 



1 



F^,„{v)dy'', 



(G19) 



where the integration is along any path connecting points 
xq and X, which belongs to the domain [/, and xo is the 
point where yl^(a;o) = 0. 
The action has the form: 



S[ti,t2] ^ - dx- dy'^F^^y)- dtH{x{t)). 



(G20) 



With respect to this structure, the connections k are con- 
stant, i.e. every function f(x,p) on the phase space com- 
mutes with them: 



f^l DB 



0. 



(G15) 



On this way, one can consistently work with the reduced 
phase space F* :ee T/{k = 0} = {(a;,p = A{x))\yx S 
M^}. The Poisson bracket on this space is defined as 



{f,9}pB on r* := {f,g}DB on rU=o ■ 



(G16) 



Since the coordinates x'^ do not commute with respect 
to the Dirac bracket 



{x^ , X }pip — , 



(G17) 



on the space F, their Poisson bracket on F* are non- 
vanishing, too: 



{x^^x''} 



PB 



(G18) 



From non-singularity of the tensor J one can conclude 
that the half of coordinates x can be used as the real 
coordinates on the F*, and rest of them are the corre- 
sponding conjugated momenta. 



The first term is the weighted surface integral over the 
surface spanned by points xq, xi and X2. The second 
term is the line term, i.e. it lives only on the line that 
connects the points xi and X2- 
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